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The charge on an atom at a metallic surface in an electric 
field is defined as the field-derivative of the force on the atom, 
and this is consistent with definitions of effective charge and 
screening charge. This charge can be found from the shift in 
the potential outside the surface when the atoms are moved. 
This is used to study forces and screening on surface atoms of 
Ag(00f ) c(2 x 2) - Xe as a function of external field. It is found 
that at low positive (outward) fields, the Xe with a negative 
effective charge of —0.093 |e| is pushed into the surface. At 
a field of 2.3 V A -1 the charge changes sign, and for fields 
greater than 4.1 V A -1 the Xe experiences an outward force. 
Field desorption and the Eigler switch are discussed in terms 
of these results. 
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The force on surface and adsorbate atoms in an ap- 
plied electric field is important for understanding field 
evaporation , field- induced reconstructions E L and the 
movement of adsorbate atoms by STM tips Qi). In the 
limit of low field, the force normal to the surface on atoms 
of type i is given by: 



(1) 



where q* is the effective charge on atom i J|-|| . £ is the 
external electric field above the surface, far away enough 
from the atoms so that it is uniform - this formula builds 
in all the local field effects and screening at the surface. 
The effective charge is also responsible for the change in 
work function <p when atoms of type i are displaced [Mb] : 
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- here Zi is the displacement of the atoms perpendicular 
to the surface, and Ai is the area per atom i. The surface 
effective charge is already familiar to us from the theory 
of the interaction of probes such as EELS and IR with 
surface vibrations jfj], and it is the same as the Born 
effective charge in polar solids in which there is much 
current interest |{L(J. In this letter we shall generalize 
these ideas to calculate the force on an atom at a metallic 
surface in an electric field of arbitrary strength. This will 



allow us to determine not only the effective charge, but 
also to assign the screening charge to individual atoms. 

We use an extension of the classical argument for find- 
ing the force on the plate of a parallel plate capacitor 
- one plate consists of the surface under consideration, 
and the other is an arbitrary electrode. The potential of 
the electrode is initially V relative to our surface, which 
is maintained at constant potential, and the two are dis- 
connected. On taking charge dq from the electrode to 
the surface the energy change of the system is dU, with: 
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If we now move the atoms of type i by dzi, the potential 
across the capacitor changes (figure 1), in just the same 
way that the work function changes in the zero-field case. 
From (3) we have: 
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But —dll/dzi is the total force on the atoms of type i, 
so (4) becomes: 
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where Fi is the force on each atom of type i, and there are 
Ni such atoms on the surface. Replacing dq by the change 
in the electric field d£ between the plates we obtain a 
relationship between dFi/d£ and dV/dzf. 
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We define dFi/d£ as the charge Qi on atoms i, and we 
can then determine it from the right-hand side of (6). 
(The difference in sign between (6) and (2) arises because 
V is the electrostatic potential, whereas 4> is the electron 
potential energy.) The right-hand side of this expression 
is easy to evaluate in an electronic structure calculation 
in the presence of an applied field - we simply have to 
move the atoms of type i and see how much the potential 
in the vacuum shifts. 
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Qi is also related to the shift of the centre of gravity of 
the screening charge with atomic displacement. The cen- 
tre of gravity z$ is the electrostatic origin of the electric 
field 0@, which means that: 
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But £ is related to the total screening charge Q per sur- 
face unit cell area A by: 



£ = 
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so from (6): 
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Now Ai I A = 1 1 Mi , where A/i is the number of atoms of 
type i per unit cell, so the charge associated with these 
atoms is given by: 
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This is an obvious way to divide the charge on the surface 
between the different atoms, and we conclude that (6) is 
a natural and unambiguous definition of Qi. 

By moving all the surface atoms through dz we simply 
shift the total screening charge, hence: 



E 



dz 
dz,. 



= 1. 



So we have the sum-rule: 



^2KQi = Q, 



(ii) 
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and the charge on all the atoms adds up to the total 
screening charge. If we separate Qi into the zero-field 
effective charge q* plus a term linear in the external field: 



Qi 



ai£, 



(13) 



we see from (6) that the force on atoms of type i in the 
field is given by: 



-a % £ 
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- the factor of 1/2 is familiar to us from the force on a ca- 
pacitor plate. In these expressions for the force, we have 
implicitly assumed that the atoms have not been allowed 
to relax from their zero-field positions; what actually 
happens is that the atoms move in the field so that this 
electrostatic force is balanced by the interatomic forces, 
in this way distributing the force as a stress throughout 
the system. 



To show how these ideas can give useful informa- 
tion about bonding and screening, we consider Ag(001) 
c(2 x 2) - Xe, with the Xe atoms adsorbed in atop sites. 
We have performed self-consistent electronic structure 
calculations for this system, using the embedding method 
]l3[ to include the semi-infinite substrate and treating 
the top two layers of atoms explicitly. We calculate the 
charges Qi for a range of applied electric fields from (6), 
by moving the atoms through small displacements (typi- 
cally 0.02 a.u.) and calculating the change in electrostatic 
potential outside the surface. Our results are shown in 
table 1, and first we note that for all the fields the sum 
of charges Qi is very close to the total screening charge 
Q calculated from (8) - the sum-rule ( |l2| ) is well satis- 
fied. The slight discrepancies are due to the fact that the 
screening is confined to the top two layers of atoms in 
our calculation. 

In the limit of zero field, we find that Xe has a small 
negative effective charge, g Xe = —0.093 |e|, which is 
largely counterbalanced by a positive effective charge on 
the underlying Ag atoms. The sign of the effective charge 
agrees with field emission experiments, which show that 
a positive field (in our convention, one which corresponds 
to positive screening charge at the surface) pushes inert 
gas atoms towards the surface |lj] . 

The negative effective charge on the Xe can be under- 
stand in terms of the dipole moment of adsorbed rare gas 
atoms and its variation with distance JlJ] . Adsorbed Xe 
has a dipole moment pointing out of the surface (positive 
charge outside, negative charge inside), as is shown by the 
decrease in work function compared with the clean sur- 
face |l6j ]. The physical origin of the positive dipole mo- 
ment lies in the fact that the vacuum side is less attractive 
to electrons on the Xe atoms than the metal side, where 
the exchange-correlation potential is more attractive jl7) . 
Now a Van der Waals treatment gives a dipole moment 
varying as (zx e — zvdw)~ 4 fllql , where zvdw is the po- 
sition of the Van der Waals plane, so as zxe increases 
and the Xe atoms are moved away from the surface, the 
dipole moment decreases. The work function increases, 
and this corresponds to a negative effective charge on the 
Xe atoms. This approach gives: 
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where A(j) is the change in work-function due to Xe ad- 
sorption. Taking (zxe ~ z vdw) = 4.2 a.u., Acf) = —0.46 
eV, and Ai = 60.7 a.u., we then obtain an effective 
charge of —0.08 |e|, in excellent agreement with our first- 
principles calculation. However, this is not the whole 
story, as a tight-binding calculation of Xe on AI shows 
some chemisorption, with real charge transfer of about 
-0.1 |e| to the Xe fl|]. Whatever the origin of the ef- 
fective charge of —0.093 |e| on the Xe, it is this charge 
which determines the force in the low field limit. 
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As the electric field is turned on, the charge on the 
Xe atoms becomes slightly less negative, whereas there 
is a much bigger increase in positive charge on the Ag 
atoms, particularly Ag(l) - the Ag atom directly under- 
neath the Xe. These results mean that the screening 
charge is mostly situated on the surface Ag atoms, in 
agreement with our expectations. What is surprising is 
that Ag(l) is responding twice as much to the electric 
field as Ag(2), and there are several possible explana- 
tions. First, the Xe lowers the Ag work-function, and if 
this is considered as a local effect we would expect that 
the valence wavefunctions on the Ag atoms under the Xe 
should penetrate deeper into the vacuum J2(J and hence 
respond more strongly to external fields. Secondly, it 
may be a local field effect with enhancement of the field 
under the Xe atoms. 

The results for the charge on the Xe shown in table 1, 
extended to larger electric fields up to £ = 0.06 a.u., can 
be fitted by the formula: 



-0.097+ 1. 1655 + 23. U7£ 2 (a.u. 
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At a field of 0.044 a.u. (2.3 V A" 1 ), the field-induced 
charge cancels the effective charge, and the charge on 
the Xe atoms changes sign. The force on the Xe is given 
by integrating ([if]) with respect to £ (6): 



-0.0975 + 0.583£ 2 + 7.715£ 3 



(a.u. 
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In a positive electric field, the Xe atoms feels a force to- 
wards the surface at low fields. The force towards the 
surface initially increases linearly with increasing field, 
but then decreases, going through zero at £ = 0.08 a.u. 
(4.1 V A -1 ) and directed away from the surface for larger 
fields. We emphasise that this is the total force on the Xe 
atoms maintained at the original position, and it takes 
into account all the field-induced bonding and charging 
plf , and local field enhancement effects. The change in 
sign of the force at large fields agrees with cluster calcu- 
lations of adsorbed rare gas atoms by Nath et al. [E2|. 
The rapid variation of fxe a t large £ means that a field 
somewhat greater than 4.1 V A -1 will overcome the Van 
der Waals and other bonding forces and remove Xe atoms 
from the Ag(001) surface. In fact this certainly overesti- 
mates the field needed for field desorption - it does not 
take account of the non-adiabatic transition to ionized 
Xe. 

The fields necessary for field desorption are much 
greater than those involved in the manipulation of sur- 
face Xe atoms by the STM. In the Eigler switch p3[ , in 
which a positively biased tip pulls a Xe atom from the 
surface and the reverse bias switches it back, typical pa- 
rameters are a tip bias of +0.8 V with the tip about 4 A 
above the surface. This corresponds to an electric field £ 
of —0.004 a.u. (in our sign convention), so from (16) we 
see that <2xe 1S essentially the effective charge. This field 
pulls the Xe towards the tip with a force of 4 x 10 -4 a.u., 



20 meV A -1 , much smaller than the force needed to pull 
a Xe atom directly off the surface (from potential energy 
curves this is 100-150 meV A^ 1 As shown by de 

Andres et al. Q and Walkup et al. ]15|, the importance 
of the force on the effective charge is that it produces a 
bias-dependent shift of the potential well of the Xe ad- 
sorbed on the tip relative to the well for adsorption on 
the surface, facilitating the Xe atom jumps. 

We thank H.J. Kreuzer for helpful correspondence, and 
R.G Forbes for useful discussions. 
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FIG. 1. Shift in electron potential energy outside the sur- 
face in an electric field of +0.02 a.u., on displacing Xe atoms 
by 0.04 a.u. into the vacuum. The undisplaced Xe atoms are 
at +2.754 a.u. 





E = 


+0.005 


+0.01 


+0.02 


Xe 


-0.093 


-0.089 


-0.084 


-0.069 


Ag(l) 


+0.086 


+0.099 


+0.115 


+0.136 


Ag(2) 


+0.002 


+0.008 


+0.016 


+0.031 


Q 





+0.024 


+0.048 


+0.097 



TABLE I. Charge on atoms at Ag(001) c(2 x 2) - Xe, for 
different electric fields. The Xe atoms are atop Ag(l). The 
last row shows the total screening charge per surface unit cell 
calculated from (8). Charges are in units of |e|, and fields in 
atomic units; a positive field is directed out of the surface. 
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